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Abstract. - We consider the damped hyperbolic equation 

euu + u t = u xx + F(u) , xGR, t > , wGR, 

where e is a positive, not necessarily small parameter. We assume that F(0) = 
F(l) — and that F is concave on the interval [0,1]. Under these assumptions, 
our equation has a continuous family of monotone propagating fronts (or travelling 
waves) indexed by the speed parameter c > c*. Using energy estimates, we first 
show that the travelling waves are locally stable with respect to perturbations in a 
weighted Sobolev space. Then, under additional assumptions on the non-linearity, 
we obtain global stability results using a suitable version of the hyperbolic Maximum 
Principle. Finally, in the critical case c = c*, we use self-similar variables to compute 
the exact asymptotic behavior of the perturbations as t — > +00. In particular, setting 
e = 0, we recover several stability results for the travelling waves of the corresponding 
parabolic equation. 
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1. Introduction 

Mathematical models for spreading and interacting particles or individuals are very 
common in chemistry and biology, especially in genetics and population dynamics. 
If the spatial spread of the particles is described by Brownian motion, these models 
usually take the form of reaction-diffusion equations or systems for the population 
densities [10,27]. Depending on the precise form of the interaction, such systems ex- 
hibit interesting solutions like propagating fronts or travelling waves, whose existence 
and stability properties have attracted a lot of attention in recent years [32]. 

It can be argued, however, that diffusion is not a realistic model of spatial spread 
for short times, because particles performing Brownian motion can move with arbi- 
trarily high speed, and the directions of motion at successive times are uncorrelated. 
This drawback can be eliminated by replacing the diffusion process with a velocity 
jump process which is more satisfactory for short times and has equivalent long-time 
properties [17,22,31]. In one space dimension, this procedure leads to damped wave 
equations instead of reaction-diffusion systems [6,19,20,33]. Under the same assump- 
tions as in the parabolic case, the damped hyperbolic equations also have travelling 
wave solutions [18] with analogous stability properties [14,16]. The aim of this paper 
is to review some of these stability results in the simplest case of a scalar equation 
with a non-linearity of "monostable" type. 

We thus consider the damped hyperbolic equation 

eu tt + u t = u yy + F(u) , (1.1) 

where y E R, t E R + and £ is a positive, not necessarily small parameter. We 
assume that the non-linearity F E C 2 (R, R) has the following properties: 

(HI) F(0) = F(l) = , F'(0) > , F'(l) < , F"(u) < for u E (0, 1] . 

In particular, u = 1 is a stable equilibrium point of Eq.(l.l), and u = is unstable. 
For simplicity, we assume F being concave on [0,1], but this condition is more 
restrictive than what we really need and could be relaxed in several ways. A typical 
non-linearity satisfying (HI) is F{u) = u — u m , with m E N, m > 2. 

Under the assumptions (HI), Eq.(l.l) has monotone travelling wave solutions 
(also called propagating fronts) connecting the equilibrium states u = 1 and u = 0. 
Indeed, setting 

u(y,t) = h(y/l + ec 2 y-ct) = h{x) , (1.2) 
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where c > 0, we obtain for the function h the ordinary differential equation 

h"(x) + ch'(x) + F(h(x)) = , xeK. (1.3) 

As is well-known [1,26], Eq.(1.3) has a solution satisfying h'{x) < for all x G R, 
h(-oo) = 1, h(+oo) = if and only if c > c* = 2 A /F 7 (0). This solution is unique up 
to a translation in the variable x. Therefore, for all e > 0, Eq.(l.l) has a continuous 
family of monotone travelling waves indexed by the speed parameter c > c*. Note 
that the actual speed of the wave is not c, but c/ Vl+£c 2 , a quantity which is bounded 
by 1 / ^fe for all c > c* . 

In the limit £ — > 0, Eq.(l.l) reduces to the semilinear parabolic equation u t = 
u yy + F(u), which has been intensively studied since the pioneering work of Fisher 
[11] and Kolmogorov, Petrovskii and Piskunov [26]. Using the parabolic Maximum 
Principle and probabilistic techniques, the long-time behavior of a large class of 
solutions has been explicitly determined [1,2]. In a more general context, a local 
stability analysis of the travelling waves using functional-analytic methods has been 
initiated by Sattinger [30] and extended by many authors [7,23,25]. In particular, 
the decay rate in time of the perturbations in the critical case c = c* has been 
computed using renormalization techniques [3,13]. Since the equilibrium state u = 
ahead of the front is linearly unstable, all these stability properties are restricted to 
perturbations which decay to zero at least as fast as h itself as x — > +oo. Following 
[14,16], the aim of this paper is to show how these results can be extended to the 
hyperbolic case e > 0. 

To investigate the stability of the travelling wave (1.2) as a solution of (1.1), we 
go to a moving frame using the change of variables 

u(y,t) = v(\/l + ec 2 y - ct,t) =v(x,t) , 

where x = VT+ £c 2 y — ct. The equation for v is 

£v tt + v t - 2ecv xt = v xx + cv x + F(v) , (1.4) 

and, by construction, h(x) is a stationary solution of (1.4). Setting v(x,t) = h(x) + 
w(x,t), we obtain for the perturbation w the equation 

ewu + w t — 2ecw xt = w xx + cw x + F'(h)w + Af(h, w)w 2 , (1.5) 
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where 

Af(h,w) = [ (1- a)F"(h + aw)da . 
Jo 

Rewriting in the usual way the second order equation (1.5) as a first order system 
for the pair (w,w t ), we shall study the stability of the origin (w,w t ) = (0,0) in an 
weighted Sobolev space Z which we now describe. 

For k G N, we denote by H fc = H fc (R) the usual (real) Sobolev space of order k 
over R, with H°(R) = L 2 (R). Following [30], we introduce for s > the weight func- 
tion p s (x) = 1 + e sx and, for k G N, we denote by the Hilbert space H fc (R, p 2 dx) 
defined by the norm 




Setting L 2 = H°, we define the product space Z — Hj, x L 2 equipped with the norm 

11(^1,^2)111 = IKHli + {Mil* ■ (i-7) 

Finally, it will be convenient to denote by Z s the space Z equipped with the s- 
dependent norm 

\\(w 1 ,w 2 ) || % e = INIlli +^11^211^2 . (1.8) 

The perturbation space Z clearly depends on the choice of s > and becomes 
smaller when s increases. It is thus natural to look for the the smallest value of s > 
for which the origin in (1.5) is linearly stable in Z. This is most conveniently done 
by setting w(x,t) = e~ sx uj(x,t) and studying the equation for u>, namely: 

eutt + (1 + 2ecs)uj t - 2ecu xt = 

(1.9) 

w« + (c - 2s)u x + (F'(h) -cs + s 2 )uj + J\f(h, e~ sx u)e~ sx u 2 . 

A straightforward computation in the Fourier variables shows that the origin in (1.9) 
is linearly stable in H 1 x L 2 only if -F"(0) — cs+s 2 < 0. In fact, this condition can easily 
be inferred from the coefficient of uj in (1.9). Therefore, the largest perturbation 
space Z for which we can expect stability of the front h is obtained by choosing 



s = -(c- \/c 2 - c 2 ) 



(1.10) 
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Note that this value corresponds to the decay rate of h as x — > +00, since h(x) ~ e _sx 
if c > c* and h(x) ~ xe _sx if c = c* [1]. Thus, if (1.10) holds, the translations 
h(x+xo) — h(x) do not belong to the perturbation space H*. In view of the translation 
invariance of Eq.(l.l), this is of course necessary to obtain an asymptotic stability 
result. In the sequel, we always assume that the condition (1.10) holds. 

In this paper, we present three different results which show that the travelling 
wave h is stable with respect to perturbations in Z e or in a subspace of it. In 
Section 2, we give a local stability result valid for all c > c* and all £ > 0. Using 
appropriate energy functionals, we show that, if (io(0), m^(0)) is sufficiently small 
in Z e , then the solution (w(t),wt(t)) of (1.5) stays in a neighborhood of the origin 
in Z e and converges to zero as t — > +00 in a slightly weaker norm. In the critical 
case c = c*, we also give an estimate of the convergence rate. These results have 
been proved in [14] in the particular case where F(u) = u — u 2 . Their proofs can be 
adapted, with minor changes, to cover the general case of a non-linearity F satisfying 
(HI). 

Section 3 is devoted to global stability results. We first recall the Maximum 
Principle for hyperbolic equations [29] in a version adapted to our problem. Then, 
under the additional assumption that F'{u) be strictly negative for u > 1, we show 
that the travelling wave h is stable with respect to "large" perturbations in Z £ , 
provided some positivity conditions are satisfied. Furthermore, if 1 + 4eF'(l) > 
and F"(u) < for u > 0, we obtain linear upper and lower bounds for the solutions 
of Eq.(1.5), as well as a decay rate in time for the quantity \\p s w(t) • The proofs 
of these results can also be found in [14] if F{u) = u — u 2 . 

In Section 4, we restrict our analysis to the critical case c = c*, and we study 
the long-time behavior of the solutions (w,w t ) of (1.5) in a slightly smaller function 
space. In particular, we show that 

w(x, t) = £f h'{x)<p* {^j^j + o(t~ 3 / 2 ) , f - +00 , 

where a G R and (p* : R — > R is a universal profile. In the parabolic case e = 0, this 
asymptotic expansion has been obtained by Gallay [13] using the renormalization 
group method [5] combined with resolvent estimates. We follow here a simpler and 
fairly different approach based on self-similar variables and energy estimates only. 
We refer to [16] for the detailed proof. 
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To conclude this introduction, we would like to point out the striking similarity 
between the stability results presented here and the corresponding statements in the 
parabolic case. This is an illustration of the more general fact that the long-time 
behavior of solutions to damped hyperbolic equations such as (1.1) is essentially 
parabolic [15]. A similar phenomenon can be observed in the context of hyperbolic 
conservation laws with damping [21,28]. 

2. Local Stability of the Travelling Waves 

In this section, we show that the travelling wave h is stable with respect to sufficiently 
small perturbations in the space Z £ . Our main result is: 

Theorem 2.1. Assume that (HI) holds, and let Eq > 0, c > c*. Then there 
exist constants S > and K > 1 such that, for all < e < e , the following 
result holds: for all ((po, (fi) G Z £ such that \\(<po, (pi)\\z e < there exists a unique 
solution (w,w t ) G C°([0, oo), Z £ ) of (1.5) with initial data (w(0),w t (0)) = ((po,<pi). 
Moreover, one has 

\\(w(t),w t (t))\\ Ze < K \\( Vo , Vl )\\ Ze , t>0, (2.1) 

and 

lim (\\w(t)\\ H , + \\(p s w(t)) x \\ L 2 + \\p s w t (t)\\ L2 ) = , (2.2) 

t — >+oo 

where p s (x) = 1 + e sx and s is given by (1.10). 

Sketch of the proof. We follow the lines of the proof of Theorem 1.1 in [14]. Let 

Ni(h) = [ F"(ah)da, N 2 (h,w) = )- f (l-a) 2 F"(h + aw) da . 
Jo 2 J 

Given < e < Eq and c > c*, we assume that (w, w t ) G C°([0, T], Z e ) is a solution of 

(1.5) satisfying 

(Al) IK*)lk ^ 5 ' te[0,T], 

for some (sufficiently small) 5 > 0. As in (1.9), we set w(x,t) = e~ sx u(x,t). To 
control the behavior of w on [0, T], we introduce two families of energy functionals: 

Eo(t) = \ ! {eu? t +uj 2 x - w 2 (/iA/"i(/i) + 2wAf 2 (h, w))) dx , 
2 Jr 

E 1 (t) = / ((i + ecs)uj 2 + eujut) dx , 
E 2 (t) = a E (t) + E 1 (t) , 
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where ao = max(2e, l/(2c 2 )), and 

£ (t) = - I {ew 2 t +w 2 x + s 2 w 2 - w 2 (hAf!(h) + 2wAf 2 (h, w))) dx , 
= / ((| ~~ £cs)ty 2 + ewwt) dx , 

./R 

£ 2 (*) = a 1 ^o(t)+^i(t)+a2(^o(t)^ 2 (t)) 1/2 , 

where «i = max(2£, #/(2c 2 )) and 9, a 2 are positive constants. Note that the func- 
tional Ei, Si control the behavior of the perturbation w "ahead" and "behind" the 
front respectively. 

Arguing as in [14, Section 2], one can show that, if S, 9 and 1 are sufficiently 
small, then the functions Eq, E 2 , £ 2 are non-negative and satisfy the differential 
inequalities 

^(t) < \{c 2 -cl)E«{t), ^(t) + E (t)<0, 

(2.3) 



d£ 2 



(t) + a 3 £ 2 (t) < C 1 (E Q (t)E 2 (t)) 1 / 2 , t E [0, T] , 



dt 

for some a 3 , C\ > 0. In addition, there exists a constant C 2 > 1 such that 

C^\\(w,w t )\\% e < E 2 (t)+£ 2 (t) < C 2 \\(w,vH)\\%,(l + nWt)h-)) , (2-4) 
where the function \1> : R + — > R + is defined by 

ty(K) = sup \F'(u)\ . (2.5) 

0<u<l+if 

Combining the estimates Eqs.(2.3), (2.4), we obtain the bound 

||(u;(f),u; t (f))|k < C ||(«;(0),^(0))|| Z£ (l + vl>(||«;(0)|| L ^)) 1 / 2 , t e [0, T] , (2.6) 

where Co is a positive constant depending only on £0, c and F. Since the Cauchy 
problem for Eq.(1.5) in Z e is locally well-posed [14], this proves global existence of 
the solution (w, w t ) provided the right-hand side of (2.6) is smaller than the quantity 
5 appearing in (Al). Then, the differential inequalities (2.3) imply that E (t) and 
£ 2 {t) converge to zero as t — > +00, thus proving (2.2). □ 
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Remarks. 

1) In the critical case c = c*, it follows from (2.3) that Eq is non-increasing in time, 
and that tE (t) + t 1 / 2 £ 2 (t) — > as t — > +oo. In particular, we have 

lim (t 1/4 (||w(t)|| H i + ||t^(t)||i^) + t 1 / 2 (||a; a (t)|| La + ||^(t)||^)) =0. (2.7) 

This decay rate is not optimal, see Section 4 below. 

2) Combining (2.1), (2.2), (2.7), we obtain 

lim ||p s «;(£)||l°° =0 if c > c* , and lim t 1 / 4 ||p s w(t) =0 if c = c* . 

(2.8) 

3) All the estimates in the proof of Theorem 2.1 are uniform in s for s G (0, £o]- m 
particular, taking the limit e — * everywhere, we obtain a proof of the local stability 
of the travelling wave h for the corresponding parabolic equation. 

4) If c = c*, the stability condition F'(0) — cs + s 2 = (s — c*/2) 2 < implies s = c*/2, 
hence (1.10) is the only possibility. On the other hand, if c > c* and s is chosen 
so that F'(0) — cs + s 2 < (in contrast to (1.10)), one can show using the same 
spectral estimates as in the parabolic case that the origin in (1.5) is exponentially 
stable in Z £ . The fastest decay rate in time for the perturbations is obtained if we 
set s = s(e), where 

c /l + 4eF'(0) 
S{£) = 2 V 1 + ec 2 • 



3. Global Stability of the Travelling Waves 

Throughout this section, we assume that the non-linearity F satisfies (HI) and 
(H2) F'(u) < -/j < , u > 1 , 

for some fi > 0. Under this additional assumption, it is known in the parabolic 
case that the front h is stable with respect to large perturbations in satisfying 
a positivity condition. This property is a consequence of the classical Maximum 
Principle for parabolic equations. In this section, we show a similar global stability 
result for e > using a hyperbolic Maximum Principle. 
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Our starting point is the following crucial observation. Let < e < Eo, c > c*, 
d G (0,1], and assume that (w,wt) G C°([0,T],Z e ) is a solution of (1.5) satisfying, 
instead of (Al), 

(A2) w(x,t) > -(1 - d)h(x) , xGR, £g[0,T]. 

In other words, we assume that v{x,t) = h(x) + w(x,t) > dh(x). Then it can be 
shown (see [14] in the case F(u) = u — u 2 ) that the a priori estimates (2.3), (2.4), 
(2.6) still hold for the solution (w,w t ), with constants Co,Ci,C2 depending only 
on e , c, d and F. Thus, we can remove the smallness condition (Al) in the proof 
of Theorem 2.1 provided we are able to show that the positivity condition (A2) is 
satisfied for all times. This in turn can be obtained from the Maximum Principle 
under appropriate assumptions on the initial data. 

3.1. The Hyperbolic Maximum Principle 

Motivated by (1.4) or (1.5), we consider the hyperbolic operator L with constant 
coefficients 

Lu = u xx + 2ecu xt - euu + cu x - u t . (3.1) 
Assume that £ : R x [0, T] — > R is a continuous function such that 

e(x,t) > I, x G R , t G [0, T] , (3.2) 

where T is a positive number and £ G R satisfies 

1 + 4e£ > . (3.3) 

The following result is a consequence of the Maximum Principle given by Protter 
and Weinberger (see [29, Chapter 4, Theorem 1] or [14, Appendix A, Theorem A.l]). 

Theorem 3.1. Let e > 0, c > 0, and assume that the conditions (3.2) and (3.3) 
are satisfied. If (u(x, t), u t (x, t)) belongs to C°([0, T],Hl oc (R) x Lf oc (R)), with u xx + 
2ecu xt - eu tt in Lf oc (R x (0, T)), and if 

(L + £(x,t))u(x,t) > , a.e.(x,t) G R x [0,T] , (3.4) 



u(x,0) < , Vi G R , (3.5) 
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eu t (x, 0) — ecu x (x, 0) H — u(x,0) < 0, a.e. in R , (3-6) 
then u(x, t) < for all (x, t) G R x [0, T]. 

As an application, we define for d G [0, 1], K > the function 

-F(u) 



A d (K) = inf 



v — u 



0<u<d, u<v < 1+K } < . (3.7) 



We have the following result: 



Proposition 3.2. Assume that (HI), (H2) hold. Let e > 0, c > c*, G [0,1] and 
Jet K be a non-negative constant such that 

l + 4eA d (K) > . (3.8) 

For some T > 0, assume that (w, w t ) G C°([0, T], Z e ) is a solution of (1.5) with initial 
data ((fiQ,ipi) satisfying 

<Po(x) > -(1 - d)h(x) , iGR, (3.9) 

e<pi(x) > ec(if' (x) + (l-d)h'(x)) - -(<p (x) + (l-d)h(x)) , a.e. in R . (3.10) 
Suppose moreover that 

w(x,t)<K, (x,t) G R x [0,T] . (3.11) 

Tien 

w(x,t) > -(1 - d)^(x) , (x,t) G R x [0,T] . (3.12) 

Proof. Without loss of generality, we may assume that F'(u) > for all u < 0. Let 
u(x,t) = dh(x) — v(x,t), where v(x,t) = h(x) + w(x,t). Since v(x,t) is a solution 
of (1.4), it is straightforward to verify that (L + l)u{x,i) = F(dh(x)) — dF(h{x)), 
where L is defined in (3.1) and 

_ F(v(x,t)) - F(dh(x)) 
{,) v(x,t)-dh(x) 

In view of (HI), one has F(dh(x)) — dF(h(x)) > for all iGR. Furthermore, since 
v(x,t) < 1 + K by (3.11), we have £(x,t) > £ = A d (K). Indeed, this inequality 
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follows immediately from (3.7) if v(x,t) > dh(x); in the converse case, we observe 
that £(x,t) > mm(F'(dh(x)),0) > A d (K). Thus (3.8) implies (3.3), and the con- 
ditions (3.9), (3.10) are nothing else as the hypotheses (3.5) and (3.6). Therefore 
Theorem 3.1 shows that u(x,t) = dh(x) —v(x,t) < for all (x,t) G R x [0, T], which 
is (3.12). □ 

Remark. Theorem 3.1 suggests the definition of a partial order in Hl oc (R) x Lj 2 oc (R) 
as follows. We say that (</?o,</?i) < (ipo,ipi) if 

¥>o(:r) < ipo(x) , x G R , 
s<pi(x) - ec<p' (x) + -(fio(x) < eipi(x) - ec^' (x) + -ipo(x) a.e. in R , 

see (3.5), (3.6). Then, if (tpo,(pi) < (V'o^V'i)) the solution of the linear hyperbolic 
equation (L + £)u(x,t) = satisfying u(x,0) = ipo(x), u t (x,0) = <pi(x) stays for all 
t G R-|- below the solution of the same equation with initial data (ipo, ipi). This order 
has the property that we can write any (</?o, <£i) G H^ oc x L^ oc as the sum of a "pos- 
itive" part {<PQ->(Pi) > an d a "negative" part (<Pq,<Pi) < 0. This decomposition 
is unique if we impose that (cpQ,(p^) = whenever (<^o,<£i) < and (<Pq,(Pi) = 
whenever (ipo, ipi) > 0. The formulas for ((Pq, tpf) are given by 

(Po(x) = sup(0,v?o(^)) , 

1 1 (3-13) 

(p+(x) = c((p$)'(x) - —iPo(x) + sup(0, {(fx - ctp' + —tp )(x)) , 

and 

(Po(x) = mf(0,(p (x)) , 

1 _ 1 (3-14) 

ip x (x) = c(tp Q )'(x) - — v? (x) + inf(0, (ipx - ap' + — (p )(x)) . 

Remark that, if (</?o ; <£i) £ -^e, then ((Pq, tpf) G Z e . Moreover, it can be verified that 
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3.2. A General Global Stability Result 

Combining the a priori estimates of Section 2 and the Maximum Principle, we are 
now able to state and prove our main stability result: 

Theorem 3.3. Assume that (HI), (H2) hold, and let e > 0, c > c*, d G (0,1]. 
There exists a constant Co > 1 such that, for all < e < Sq and all K > satisfying 

l + 4eA d (K) > , (3.15) 

the following result holds: If > is such that 

CqK*{1 + ^(if*)) 1 / 2 < K , (3.16) 

where ^> is defined in (2.5), then for any (<p , <Pi) G Z £ satisfying the inequalities (3.9), 
(3.10) and the bound \\(<fio, <fii)\\z s < K*, there exists a unique solution (w,wt) G 
C°([0, oo), Z e ) of (1.5) with initial data (<po, (pi). Moreover, one has 

\\(w(t),w t (t))\\ Ze < K , w(x,t) > -(l-d)h(x) , (3.17) 

for all x E R, t G R+, and (2.2), (2.7) hold. 

Remark. The constant C in (3.16) is the same as in (2.6). 

Sketch of the proof. By Proposition 3.2, the solution (w, Wt) with initial data 
(</?o, ^l) satisfies w(x, t) > — (1 — d)h(x) as long as w(x, t) < K. On the other hand, 
in view of (3.16), the a priori estimate (2.6) implies that w(x, t) < \\(w(t),wt(t))\\z £ < 
K as long as w(x, t) > —(l-d)h(x). Combining these facts and using a contradiction 
argument, we show that the solution is globally defined and satisfies (3.17) for all 
times. The differential inequalities (2.3) then imply (2.2), (2.7). □ 

Remarks. 

1) The relations (3.15), (3.16) imply that K (hence K*) can be chosen very big if 
e is sufficiently small. In this case, Theorem 3.3 shows that the travelling wave h is 
stable with respect to large perturbations, provided the positivity conditions (3.9), 
(3.10) are satisfied. Conversely, if e is large, then K (hence K*) has to be very small, 
and Theorem 3.3 reduces to a local stability result similar to Theorem 2.1. 
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2) As a simple example, consider the non-linearity F(u) = u — u 2 [14] which satisfies 
(HI), (H2). Then Ad(K) = -(d+K), and the condition (3.15) reads 1 -Ae(d + K) > 
0. Also ^(K) = 1 + 2K, hence the hypothesis \\(<fio, <fii)\\z e < K* can be replaced by 

\\{^^i)\\z £ <CK{l + K)- l /\ 

where C is a sufficiently large positive constant. 

3.3. Further Stability Results if 1 + 4eF'(l) > 

The previous results are still incomplete for at least two reasons. First, if c > c*, 
they fail to give a decay rate of the perturbations as t — > +oo, see (2.8). Next, they 
do not provide any global existence result if d = 0, that is if v(x,0) > 0. Assuming 
that the non-linearity F satisfies the stronger assumption 

(H3) F"(u) < , u > , 

we can give a partial answer to both questions when 1 +4eF'(l) > 0. Indeed, in this 
case, the Maximum Principle allows us to compare the solution w(x,t) of (1.5) with 
solutions of the linear equations 

ew u +w t - 2ecw xt = w xx + cw x + F'(h)w , (3.18) 

and 

ew u +w t - 2ecw xt = w xx + cw x + F'(0)w . (3.19) 

We denote by (w(t),w t (t)) the solution of (3.18) with initial data ((po,(pi), and by 
(w ± (t), wf(t)) the solution of (3.19) with initial data (<p^,<pf), the "positive" and 
"negative" parts of (ip ,ipi) defined in (3.13), (3.14). Following the arguments in 
[14, Section 4], we obtain our last stability result: 

Theorem 3.4. Assume that (HI ), (H3) hold and that l+4eF'(l) > 0. Given c > c*, 
d G [0, 1], there exists a constant C 3 (c) > 1 such that, for all K > satisfying 

1 + > , 

the following result holds: for any (<£o>¥>i) £ %e satisfying (3.9), (3.10) and 
inf (ll^o, Vi)lk, \M,<pt)\\z e ) < t^t, 
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there exists a unique solution (w,w t ) € C°([0, oo), Z e ) of (1.5) with initial data 
{ipo.ipi). Moreover, we have 

-(l-d)h(x) < w(x,t) < K , (3.20) 

and 

w~(x,t) < w(x,t) < w(x,t) < w + (x,t) , 
for alJi6R,lGR + . Finally, if d > and 1 - 4eF'(0) > 0, one has 

t lim t 1 / 4 (||p fl «;(t)|| L oo + |Kt)|| H i + \\w t (t)\\v) = . 

Remark. If {<Pq,<Pi) = 0, i.e. if the initial data are non-positive, then (3.20) with 
K = shows that the solution w(x,t) remains non-positive for all times. 

4. Asymptotic Expansions in the Critical Case c = c* 

In this section, we restrict ourselves to the critical case c = c*, and we consider 
perturbations of the travelling wave h in a strict subspace of Z e . Using self-similar 
variables and energy estimates, we are able to compute explicitly the long-time 
asymptotics of the perturbations as t — > +oo. In particular, we recover the results 
obtained by Gallay [13] in the parabolic case e = 0. 

Following Kirchgassner [25], we consider solutions of (1.4) of the following form 

v(x,t) = h(x) + h'(x)w(x,^£j , 

i.e. we set w(x,t) = h'(x)W(x,t/(l+ecl)). Then W satisfies the equation 

v W tt + (l-in(x))W t -2uW xt = W xx + 1 (x)W x + h'(x)W 2 U(h(x),h'(x)W) , (4.1) 
where 

e ec* , . n h"(x) _ 

(1+ecj) 2 1 + eci h'(x) 

In (4.1) and in the sequel, the second argument of the function W is simply denoted 
by t, instead of t/(l+ec%). 
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From [1] we know that the travelling wave h (with c = c*) satisfies 

f 1 + O(e^) as x -> -oo , 

t (aix + a 2 )e c * x/2 + C(x 2 e c * x ) as x — > +oo , 



where ai > 0, a 2 G R, and /5 = |(— c* + a/c 2 — 4F'(1)) > 0. Similar asymptotic 
expansions hold for the derivatives h', h" , hence 



[ 7_ + (^(e^) as x — »■ — oo , 

as x — > +oo , 



7(X) i 2/(x+x ) + ©(xe- 6 **/ 2 ) ( L2) 



where 7_ = c* + 2/3 = 2 A /F'(0) - F'(l) and x = (a 2 /a 1 - 2/c*). The hypothesis 
(HI) on F also implies that 7 ; (x) < for all x G R. 

To study the long-time behavior of the solution W of (4.1), we use the scaling 
variables or self-similar variables defined by 

e = ^, r = iog ( wo), 

for some to > 0. These variables have been widely used to study the long time 
behavior of solutions to parabolic equations, in particular to prove convergence to 
self-similar solutions [4,8,9,12,24]. In [15], it has been shown by the authors that 
these variables are also a powerful tool in the framework of damped hyperbolic 
equations. Following [15], we define the rescaled functions U and V by 

r) = e 3 ^ 2 W(^/ 2 , e T - t ) , Vfo r) = e 5 ^(£e^ 2 , e T - t ) , (4.3) 

or equivalent ly 



x 



WM = jiT^ u {vm' logit+to)] • 



(4.4) 



Then {/(£, r), V(£, r) satisfy the system 

u T -lu,- 3 -u = v, 



- i Vi - \V) + (1 - ^e^))V - 2ue~^V, = ^ 
% + e T/2 7 (^/ 2 )^ + e-/ 2 h'(^/ 2 )U 2 N(^ r) , 
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where N(£,t) = Af(h(£e T / 2 ), e- 3T / 2 h'(£e T / 2 )U). 

We now introduce function spaces for the rescaled perturbations U, V. For 
t > 0, we denote by X T , y r the Hilbert spaces of measurable functions on R defined 
by the norms 

/0 /-+oo 
e 2 ^ T/ >(e)i 2 de+ / (i+ninoi 2 ^, 
-oo JO 

m\ T = \\u\& r + \\ut\\i r . 

We denote by Z T the product space Z T = X T x X- equipped with the standard norm 

As is easily verified, if the functions (U, V) and (W, W t ) are related through (4.3) 
or (4.4), then (U(-, r), V(-, r)) G Z r if and only if the actual perturbation w(x, t) = 
h'(x)W(x,t/(l+ecl)) satisfies 

(w 2 + wl + w 2 )(x,t) dx+ / (l + x 4 )e c * x (w 2 + wl+w 2 )(x,t)dx < oo . (4.6) 

-oo J0 

Therefore, the perturbation space considered in this section is slightly smaller than 
the space Z defined in (1.7), due to the factor (1 + x 4 ) in (4.6). 

Before stating the main result of this section, we explain its content in a heuristic 
way. Taking formally the limit r — > +oo in (4.5) and using (4.2), we see that U 
satisfies the linear parabolic equation 

U T = C^U = %+ (i + ^Us + ^U if£>0, Ut = if£<0. (4.7) 

Therefore, it is reasonable to expect that the long-time behavior of the solutions 
to (4.5) will be determined by the spectral properties of the operator on R + , 
with Neumann boundary condition at £ = 0. Now, as is easily verified, this limiting 
operator is just the image under the scaling (4.4) of the radially symmetric Laplacian 
operator in three dimensions. Indeed, if U and W are related through (4.4), the 
equation U T = CooU is equivalent to W t = W xx + (2/x)W x , x > 0. This crucial 
observation allows to compute exactly the spectrum of in various function spaces, 
see [15]. For instance, in the space L 2 (R + , (l+£ 6 )d£), the spectrum of consists 
of a simple, isolated eigenvalue at A = 0, and of "continuous" spectrum filling the 
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half-plane {A G C | Re A < —1/4}. The eigenfunction corresponding to A = is the 
gaussian e - ^ / 4 . Therefore, we expect that the solution U(£, r) of (4.5) converges as 
t — > +oo to a<p*(£) for some a G R, where 

*m - 1 J 1 if ^ <0 ' 
v {0 ~7^W e/A if € > o . 

This function is normalized so that J Q °° £ 2 <£>*(£) d£ = 1. Since V = U T — — 
we also expect that V(£, r) — > where = — |</7| — |</?*. It is crucial to note 

that Eq.(4.7) is independent of e: this explains why the solutions of (4.1), hence of 
(1.1), behave for large times in a similar way to those of the corresponding parabolic 
equations. 

Our last result shows that the heuristic arguments above are indeed correct: 

Theorem 4.1. Assume that (HI) holds, and let e > 0, c = c*. There exist to > 
and do > such that, for all {Uq, Vq) G Z To satisfying \\(Uo, Vo ) 1 1 < So, the system 
(4.5) has a unique solution (U, V) G C°([r , +oo), Z T ) with (U(t ), V(t )) = (U , V ). 
In addition, there exists a* G R such that 

\\U(T)-a*<p*f XT + f e-^/ 2 \\V(a)-a*r\\ 2 y„da = 0(r 2 e^/ 2 ) , 

J To 

as t — > +oo. 



Remark. We say that (U, V) G C°([r , +oo), Z T ) is a solution of the system (4.5) if 
there exists a solution (W, W t ) G C°([0, +oo), Z ) of (4.1) such that (4.3), (4.4) hold, 
with t = e T °. 



In terms of the original variables, Theorem 4.1 implies that [16] 

'x^/l+ecj 



^( 1 + ^)^>-^'<M^) 



o(t- 7 /Ho g t) , 



as t — >■ +oo, where a = a*(l + ec 2 ) 3 / 2 , w(x,t) = h'(x)W(x,t/(l+ecl)) and W is 
given by (4.4). In the parabolic case e = 0, this result has been obtained in [13] 
using slightly different function spaces. Remarkably enough, the asymptotic profile 
if* is universal: it does not depend on the initial data, nor on the parameter e > 0, 
nor on the precise form of the non-linearity F. 
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